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INTRODUCTION 


The need to develop large solar-cell arrays for space applications has stimu- 
lated a great interest in thin-film photovoltaics because of their projected high 
power-to-weight ratio (greater than 1 kW/kg) . Specifically, it would be advantageous 
to develop cost-effective thin-film GaAs solar cells because they possess a wide 
range of desirable characteristics for space applications (ref. 1). Such cells could 
be launched and deployed by the Space Shuttle as solar arrays to be used as long-term 
power supplies for various Shuttle- and satellite-related missions. 

A promising approach for developing these solar cells is to epitaxially grow 
high-efficiency photovoltaic thin films on low-weight, inexpensive foreign substrates 
using vapor deposition techniques. In order to construct high-ef f iciency solar 
cells, it is necessary to produce well-ordered, defect-free GaAs epitaxial films. 

The degree of crystalline order achieved will be controlled by the adsorption, nucle- 
ation, and lateral-growth behavior in the first few layers deposited. This behavior 
is, in turn, controlled by such parameters as the deposition rate, surface tempera- 
ture, adatom mobility, adatom-substrate chemical reactivity, surface defect density, 
and lattice registry of the system. In order to assess the relative sensitivity of 
these parameters and to model experimental research efforts^ a solid-on-solid (SOS) 
Monte Carlo computer simulation for a one-component system is presented that utilizes, 
a potential-energy scaling technique over a 400-site array. 

Although numerous Monte Carlo models exist (refs. 2 to 5) , the advantage of this 
approach is that it allows single events to occur simultaneously at all sites of the 
array within the constraints of the surrounding potential-energy field and the ther- 
mal fluctuations associated with a given substrate temperature. This is in contrast 
to past models, which utilized rate equations to determine multiple events for a 
given atom (in one sampling interval) while freezing the motions of all other atoms 
in the array. Continuous updating of the potential energy over all sites of the 
array on a time scale sufficiently small to cover single events at each site (e.g., 
adsorption/desorption of atoms and single-atom jumps to neighboring sites) is made 
possible by the use of ordered statistics. Therefore, the simultaneous changes at 
all sites and the resulting individual atomic interactions provide the conditions 
necessary to simulate thin-film growth in a near physically real manner. 


MODEL CONCEPTS 

The model employs a 20 x 20 square array with periodic boundary conditions as 
the substrate (and reference plane) on which atoms are deposited. As atoms are ran- 
domly added to the sites they are free to desorb, migrate, diffuse into (incorporate) 
and out of the substrate, or remain localized. These changes alter the stacking 
height of each column growing on a particular site, and also change the potential 
energy at that site and at neighboring sites. Figure 1 illustrates a typical surface 
topography (terrace-ledge-kink) and represents the surface model used for the solid- 
on -so lid (SOS) method of studying crystal growth. In this approach, atoms are repre- 
sented by cubes and are constrained to grow in well-defined columns. 



Adatom 



Figure 1*- Geometry of terrace-ledge-kink SOS model. 

The physical constraints that determine the temporal behavior of every atom 
located at the surface of each column site (ifj) are based on the interaction poten- 
tial of each surface atom relative to the substrate and to its nearest neighbors. 

The interaction potential across a perfectly homogeneous surface with well depth 

and interatomic distance a^ is depicted as changing uniformly from site to site 
(fig. 2(a)). (The symbols used in this paper are defined in a list after the refer- 
ences.) A typical interaction potential across a heterogeneous surface is shown in 
figure 2(b) . The variation in potential energy at different sites reflects the many 
surface features shown in figure 1, such as vacancies, steps, kinks, and clusters of 
various sizes. In this model it is assumed that all incident atoms adsorbed on the 
surface thermally accommodate to the surface instantly. The sticking coefficient 
S for the substrate and the growing film is governed by 


S = S (1 - 0)^ (1) 

O 

where 0 is the coverage is the sticking coefficient at zero coverage and ^ 

is the order of bonding. 

Once an atom is adsorbed onto, or migrates to, a particular site (i,j), the 
potential interaction energy at that site immediately changes from that of the atom 
below to that of the newly adsorbed atom. If an atom is desorbed or migrates away, 
the potential at that site changes to represent the atom that was immediately below 
it, that is, the ”new surface atom." Since the heat of adsorption at zero coverage 
reflects the number of bonds that an adsorbed atom makes with the substrate (e.g., 
three bonds for (111) face-centered cubic (fee) structure, four bonds for (100) fee, 
and five bonds for (110) fee), only the lateral interactions of the growing first 
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(a) Variation over a uniform homogeneous surface. 



(b) Disturbances in surface potential caused by vacancies, dislocations, 

kinks, and impurities • 

Figure 2.- Variations in potential energy at surface. 


layer need to be considered individually. In addition, the changes in potential 
energy at a particular site will affect the potential energy of the nearest neigh- 
bors, second-nearest neighbors, and third-nearest neighbors in the first layer of 
growth. Updating the potential energy at a given site following a kinetic event 
therefore requires simultaneous updating of the potential energy for all surrounding 
sites. 


The changes in potential energy that occur along a lineal section of a homoge- 
neous surface upon the adsorption of a single adatom at site (i,j) are shown in 
figure 3. Note that the potential at that site increases by an amount <j) , and the 
potential at adjacent sites decreases by an incremental amount a. Ihis represents 
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Figure 3.- Changes in potential energy along lineal section of a homogeneous 
surface upon adsorption of a single adatom at site (i,j)* 


the actual physical condition of the particular adatom at site (i/j); it requires 
less energy to desorb or migrate than did the original surface atom at that site when 
it was surrounded by all its nearest neighbors. Further, the deeper potential at the 
neighboring sites reflects the increased energy necessary to desorb an atom at those 
sites. Ihis is a result of the increase in number of neighbors for those sites. 

Figure 4 shows the variation in potential energy at site (i,j) as lateral neigh- 
bors are progressively added. For each nearest neighbor, a potential-energy incre- 
ment a is subtracted, and increments of p and y are subtracted for second and 
third neighbors, respectively. Autoepitaxy requires that the initial xx>tential 
be recovered when all the surrounding-neighbor sites have been 
filled. Thus, for an adatom at site (i,j), 

U(i,j) = -U + (t) = -AH ^ (2) 

o o ads 

For a (110) fee surface, the addition of two nearest neighbors (-2a), two second- 
nearest neighbors (-2p), and four third-nearest neighbors (-4y) to the adatom 
recovers the potential of the original site. Since 


4) = 2(x + 2p + 4y 
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Adatom-surface distance, z 


Figure 4.- Variation in potential energy at site (i,j) as lateral neighbors 

are progressively added* 

then 


U(i,j) = -u + ()) - 2a - 2p - 4y = ~U 

o o ^ ' o 

The values of a, p, and y are determined by the lateral bond strength of a 
particular system* 




Energy Fluctuations and Potential Scaling 


The behavior at each site in the entire array is monitored once during every 

sampling interval Atg • By random-number generation, the fluctuation in vibrational 

(thermal) energy e is added to the interaction potential U(i,j) at each site, and 

the resulting total energy E(i,j) for the site is compared with the magnitude of 

the energy barriers for desorption, surface migration, incorporation, and diffusion 

of a substrate atom into the film* If E(i,j) of the atom exceeds one of these 

barriers but is less than another, then the atom moves accordingly* If E(i,j) is 

less than the lowest barrier, then the atom remains localized* For each sampling 

interval At , an atom in an arbitrary site (i,j) has the total energy 
s 


E(i,j) = U(i,j) + e 


(3) 


where e is determined by random number R from the Boltzmann distribution 


f ( e) - \ exp(-Xe) 


(e > 0) 


(4) 


where \ = 1/kT* The mean energy of this distribution is 


■I 


<e> = I e f(e) de = kT = -^ 


(5) 


and the cumulative energy distribution is given by 


{ 


F(e) = I f(e) de = 1 - exp(-Xe) 
'0 


( 6 ) 


A random variate e can be generated by using the inverse function associated with 
the cumulative distribution (ref* 6)* Ihus, for the uniform random number R 
between 0 and 1 and with R = F(e), the inverse function gives the random energy 


£ = -kT ln(1 - R) 


(7) 


Equation (1) then becomes 


E(i,j) = U(i,j) - kT ln(1 - R) 


( 8 ) 


Selection of a random number during each At allows a value E(i,j) to be deter- 
mined. Figure 5 represents the ground state^potentia 1-energy scale for the adsorbate 
atoms, to which E(i,j) is compared* Basically, four cases are considered* 
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Interaction 

potential, 

U 



^ Free atom E(i, j) > 0 



Migration 

^ or incorporation o > E(i, j) > U. 




Migration > U 




i Localization 

I “m " i> 



Figure 5.- Total energy of site (i^j) at a given instant compared with threshold 

levels of atomic events. 


1. If E(i,j) > 0, then the atom desorbs and disappears from the surface. 

2. If 0 > E(i,j) > then the atom migrates or is incorporated into the bulk 

and disappears from the surface. 

3. If > E(i,j) > then the atom migrates to a vacant nearest-neighbor 

position, 

4. If Uj^ > E(i,j), then the atom remains localized. 

In the first case, the desorbing atom vacates the site and is no longer moni- 
tored. The frequency of desorption v^es dependent on the number of bonds to 
other atoms in the surface and is given by (ref. 5) 

where t) is the number of bonds, (j) is the ^^ngle-bond energy, and is the 

vibrational period of surface atoms (t « 10 sec). 

o 

In the second case, if an atom has sufficient energy to incorporate then it also 
has sufficient energy to migrate. Whether it will incorporate or migrate depends on 
the relative probability of each possibility. In the case of migration, the avail- 
ability of more than one site for migration is determined by the coordination number 
associated with the available sites. The higher the coordination number of an 
available site, the deeper the potential well and the smaller the energy barrier to 
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migration. This increases the probability of migration to that site. Ihe experi- 
mental evidence to support this assumed behavior is sizeable (ref. 7). 

Migration, or surface diffusion, can be localized or nonlocalized, depending on 
the amount of thermal energy absorbed. The case of localized adatom diffusion is 
represented by (ref. 8) where a^ is the interatomic distance (A = 1/4 for 



a (100) fee structure) and is the energy of migration. (See fig. 3.) An adatom 

in this state can be characterized by three degrees of freedom, two vibrational and 
one translational, and moves over single saddle points. If the adatom has signifi- 
cantly more energy than but less than the desorption energy, then nonlocalized 

diffusion will occur. The atom is then characterized by one vibrational and two 
translational degrees of freedom. In this case, the atom may take several hops along 
potential-energy saddle point paths, or it may leap completely over the top of a 
nearest-neighbor atom. In general, the obtained from experiments probably 

represents some combination of localized and nonlocalized diffusion. 


In the case of incorporation, an atom is considered to have been removed from 
the surface in the same way that an adatom desorbs. It is assumed that the bulk 
vacancy concentration is sufficient to receive the adatoms, and thus the atom simply 
disappears from the surface site. Hie defining equation for the diffusion of atoms 
from the adlayer (adsorption layer) into the bulk is (ref. 9) 



where C|^ is the bulk concentration in the selvedge, is the deposition rate, 

and erfc is the error function complement. The amount Ac incorporated, then, is 
the total amount absorbed by the bulk per unit area, as follows: 


Ac 

s 





( 12 ) 


In the third case, the atom has only enough energy to migrate. Surface migra- 
tion is treated as desorption followed by adsorption at a neighboring site. In the 
fourth case, no change occurs at the site. 


Initially, similar determinations must be made for each adatom that adsorbs onto 
a substrate of a different species. These determinations are made by replacing 
E(i,j) with Eg(i,j) in the inequality argument presented previously. Each inter- 
action is considered separately until all substrate sites are covered. 

After the substrate is completely covered by the adlayer, diffusion of substrate 
atoms through the growing film to the surface can be approximated by randomly depos- 
iting substrate atoms onto surface vacancy sites at a rate controlled by the diffu- 
sion equations for the adlayer-substrate system. Since the growing film is so thin. 
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this can be considered diffusion through a membrane, and the total nximber of sub- 
strate atoms emerging at the surface can be calculated by (ref. 9) 


s ^sub 




(13) 


Ordered Statistics 

The residence time or stay time, of an atom on a surface is given by the Frenkel 
equation 


T 


T 

O 




(14) 


where ^ads heat of adsorption and Tq is the period of vibration for 

surface atoms ( « 10"^^ sec). (It is assumed here that = AH^^g • ) Ideally, 

then, the most physically real sampling time that corresponds to changes in vibra- 
tional energy, and therefore to changes in E(i,j), is Atg = Computer costs, of 

course, prohibit the extensive run times that would be required to sample 400 sites 
10^^ times each second. In order to circumvent this difficulty, the method of 
ordered statistics is applied (ref. 6). Essentially, most of the time -dependent 
energy variation at a particular site results in insufficient thermal fluctuation for 
atomic movement, so the atoms are localized over most of the sampling interval. This 
large time interval of atomic localization is not important to the actual thin-film 
growth because there are no atom movements; therefore, only the fraction of the 
sampling interval in which movement does occur needs to be considered. Thus, we must 
determine what fraction of the time the site energy is in excess of the minimum 
activation barrier for movement. 


Let e1 , e2, ..., en denote n random samples from the exponential distribu- 
tion in equation (4) and let ^(i)f ^(2)' ^(n) denote the ordered arrangement 

(from low to high) of the n random samples with E(i_i) < for all 

i = 2, 3, ..., n and e(j^) = max(e1, g 2, ..., en). The probability distribution 

of the largest ordered statistic e = e, . is given by 

(n) 

g(e) = n[F(e)]’^"^ f(e) (0 < e < ») (15) 


where f ( e) and F(e) are given by equations (4) and (6). 


Figure 6 presents examples of this distribution for a temperature of 500 K and 
different sampling sizes (n = 10^, 10®, and 10®). The cumulative frequency distribu- 
tion is given by 


G(e ) 
o 


= P(e < 


e ) 
o 



— — ' 

/ g(e) <3e = 

1 - exp(-\e ) 
o 


(16) 
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€, eV 


Figure 6.- Distribution function of Boltzmann-ordered statistics at 
for different values of the sampling interval. 


T = 500 K 


To generate random variables £(nl this distribution, a uniform random number 

R is generated such that J = R. This gives 


£, , = -kT ln( 1 - 
(n) 


(0 < R < 1) 


(17) 


which can be compared with equation (6) . Note that for very large n, we can approx- 


imate e, . by 
(n) 


e, . = -kT In - — In R 
(n) \ n J 


(18) 


hops is T, 


Atg = 


10 


- 5 " 


If, for example, T = 500 K, = 0.7 eV, and the average mean time between 

= 10“^^ exp(\Q^) = 10”^ sec, then in a sampling interval of 
sec, a single hop occurs. Any smaller sampling interval is unnecessary 
because no movement occurs. Any larger sampling interval (i.e., greater than 
multiples of 10”^ sec) would result in multiple events for a single atom and multiple 
events at other sites in the array. Therefore, the more the sampling interval 
exceeds the time of a single event, the less physically real the model becomes. 


After selecting a particular system for adsorption, the minimum activation bar- 
rier determines the sampling interval and therefore the number of random samples n. 


10 



The ordered statistic £(n) determined by a random-number generator and 

equation (17), Combining equations (17) and (3) gives the total energy at site (irj) 
as 


E(i,j) = U(i,j) - kT ln(l - 


(19) 


The same procedure described previously to determine atomic events is then followed. 


PROCEDURE 


The changes in potential energy over a uniform face-centered cubic surface can 
be depicted by a mnemonic mask. When the center of this mask is placed over 
site (i,j), the mask characterizes the changes in potential energy that occur both at 
the central site and at the nearest-neighbor sites upon the adsorption of an atom. 

As shown in figure 3, the adsorption of an adatom at site (i,j) drives the potential 
energy at that site positive by an amount and simultaneously drives the poten- 

tial at the nearest-neighbor sites more negative. As other adatoms arrive at 
nearest-neighbor sites, the potential at site (i,j) incrementally becomes more nega- 
tive. In the case of a (100) fee surface, there are four nearest neighbors (incre- 
mental change a) and four second-nearest neighbors (incremental change p) in the 
surface layer. As these sites are filled, the potential at site (i,j) returns to 
Thus, for the different crystal faces 


4)^ = 4a + 4p 
for (100) fee surface. 


4)^ = 6a' + 6P' 



( 20 ) 


( 21 ) 


for (111) fee surface, and 

4)^ = 2a” + 2p” + 4y 
for (110) fee surface. 



a(a ) ; 
o 



( 22 ) 


Figure 7 gives the actual atomic arrangement for the three low-index crystal 
faces for a face-centered cubic structure, and shows how the potential-energy incre- 
ments are considered in terms of the SOS model. The potential at each site is 
assessed independently by moving the mnemonic mask to the site and evaluating all the 
positions on the mask to determine the new potential. This ultimately determines 
whether an atom has been added, lost, or unchanged since the last survey. In the 
case of atoms in a crystal that are bonded by van der Waal's forces, the effects of 
second- and third-nearest neighbors are quite significant, but for metal, ionic, and 
covalent bonds, only the primary bonds need to be considered (p = y == 0). The (110) 
face, however, is an exception (y = 0) . 
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Crystal 

face 


Adatom 


o 


Lateral- 

neighbor 

atom 


Mnemonic mask 


( 111 ) 



-iS 

-P -a -a -/3 

-a (p^ -a 

-i3 -a -a -/3 

-P 


(100) 



-Q -P 
-a 0Q -a 
-13 -a -13 



-y -a -y 


-/3 

-y -a -y 


Figure 7*- Mnemonic mask used to update potentials at site (i,j) and neighboring 

sites following the motion of an atom. 


The flow chart for this simulation is presented in figure 8. Initially, the 
boundary conditions can be established by adjusting the site potentials to reflect 
desired distributions of kink sites, ledges, impurities, or other defects. Any 
degree of heterogeneity can be included. Ordinarily, an assumed impingement flux 
dictates the number of adatoms arriving at the surface during a prescribed sampling 
interval. Depending on the value of S, adatoms are added to sites by random-number 
generation, and the mnemonic mask is then used to update the potential energy at each 
site as well as that at each of its surrounding neighbors. Every site in the array 
is surveyed, and a random thermal energy e(R) is generated and ultimately used to 
determine the total energy at site (i,j). (See eq. (19).) For each site it is then 
determined whether desorption, migration, incorporation, or localization has occurred, 
and the atom is moved accordingly. Ultimately, all the sites are independently and 
randomly sampled. This interaction procedure is followed until the desired growth 
time has expired. At this point, crystal growth measures such as vertical growth 
rate, surface diffusion coefficient, desorption rate, incorporation rate, nucleation 
rate, and cluster size distribution can be calculated. Further, graphic displays of 
surface migration paths and nucleation island growth and coalescence can be plotted, 
along with isometric views of the growth. 
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Figure 8.- Flow chart for Monte Carlo SOS computer simulation of thin-film growth. 
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RESULTS AND DISCUSSION 


Five computer experiments were conducted on a (1 00) fee surface to assess the 
physical behavior of the model* The experiments and their parametric variations are 
summarized in table I* In order to model a thin-film system, several a priori 


TABLE I.- MODEL EXPERIMENTS AND PARAMETERS 


Experiment 

-2 -1 IS 

R^, cm -sec X 10 

T, K 

t, sec 

> 

0 ) 

1 - Surface diffusion of 

0 

300, 400, 500, 

2 

0.7 

single adatom 


600, 700 



2 - Clustering of nine 

0 

400, 500, 600 

2 

.7 

dispersed atoms 





3 - Thin-film growth with 

.5 

550 

2 

.7 

a single defect site 





4 - Thin-film growth with 

in 

• 

in 

o 

. 

300, 400, 500, 

2 , 15 

.7 

variable deposition rate 

1, 2.5 

600, 700 



and substrate temperature 





5 - Annealing of thin-film 

0 

500, 600, 700, 

6 

.7 

growth 


800 




parameters must be specified. Since there are insufficient data on any particular 
semiconductor -metal thin-film system, nominal parameters were chosen for the experi- 
ments. Ihese parameters are presented in table II. 


TABLE II.- NOMINAL PARAMETERS 


<^Hads)0=O' 



D_^ (self -diffusion into growing film), cm -sec .... 

D (diffusion into substrate), cm^-sec“^ 

D (substrate atom diffusion through film), cm -sec 

2d' 

s 


1.7 

3.87 

7.8 exp(-2.98\) 
0.44 exp(-2.09\) 
0.1 3 exp(-1 .08X) 

0.7 

1 


The growth kinetics are quite sensitive to the magnitudes of these parameters. 
For example, the selected values for ^0 high, and for moderate 

to low temperatures this will yield capture-controlled kinetics. For the same tem- 
perature range, reducing ^^^ads^ 0 -O one-half will yield desorption-controlled 
kinetics. An excellent discussion of this argument can be found in reference 2. 


14 







Experiment 1 - Surface Diffusion of Single Adatoms 


The random migration of atoms over a perfectly uniform surface (fig* 2(a)) is 
essentially a random walk to nearest-neighbor sites for the (100) fee crystal face, 
which has four primary sites. The number of migrations as a function of time is 
linear and varies exponentially with temperature, A plot of the calculated diffusion 
coefficient as a function of inverse temperature is given in figure 9, The linear 
behavior of the diffusion coefficient follows equation (10), and the slope of the 
Arrhenius plot yields the activation barrier = 0,7 eV (the input condition). 

This provides a self-consistent check on the migration physics of the computer 
model. The average number of migrations for each temperature is given adjacent to 
each data point in figure 9, In the case of a nonuniform surface, both mean dif- 
fusion length and migration frequency would be substantially reduced because of the 
lower probability of escape from kink sites, steps, and other defects. 


Experiment 2 - Clustering of Nine Dispersed Atoms 

This experiment demonstrated the clustering of atoms to the most stable configu- 
ration, Each atom should ultimately maximize its number of nearest neighbors, thus 
making U(i,j) more negative and therefore more stable. Ideally, the atoms should 
interact to form a 3 x 3 array. The dispersed atoms individually random walk and 
collide to form dimers, trimers, and finally a nine-atom cluster. Figure 10 presents 
a graphic sequence of the atom motion as a function of time for a substrate at 600 K, 
In figure 10(b), many migrations have already occurred and the atoms have formed a 
monomer, two dimers, and a tetramer. Figure 10(c) shows the effective rotation of 
the tetramer and the clustering of the two dimers and the monomer. In figures 10(c), 
10(d), and 10(e), adequate thermal energy was absorbed to break away a single atom 
from the tetramer, trimer, and dimer, respectively. Eventually each monomer was 
captured by the growing cluster. Finally, in figure 10(f), the cluster has reached 
its most stable state. Cluster growth is believed to occur either by this method of 
monomer capture or by actual en masse motion of one of the clusters, subsequent col- 
lision with another cluster, and reorientation of the atoms to registry with this 
cluster (ref. 2, chapter 3), In the formulation of this model, no consideration was 
given to the motion of whole clusters, but, as these experiments show, whole clusters 
can move by individual atom motion at the periphery of the cluster, which results in 
the net motion of its center of mass. 


Experiment 3 - Thin-Film Growth With a Single Defect Site 

A generalized point defect is modeled by assuming that the potentials for the 
"trap" and its nearest-neighbor sites vary in the lateral and vertical directions and 
are dependent on layer height, (Potentials return to normal after two layers of 
growth.) Although it is understood that will vary in the vicinity of defects, 

there is no simple means of determining its value near a particular defect. It will 
therefore be assumed that = 10“^^ sec. 

Figure 11 shows the prescribed variation in potential energy at the trap, which 
is positioned at the center of the 20 x 20 array. A graphic display of the growth 
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Figure 9.- Arrhenius plot of single-adatom diffusion coefficient over a 
uniform surface as a function of reciprocal temperature. 





(a) t = 0. 


(d) t = 0.5 sec. 



(b) t = 0.1 sec. 


(e) t = 0.7 sec. 



(c) t = 0.3 sec. 


(f) t = 1 .0 sec. 


Figure 10.- Clustering of dispersed atoms on a uniform surface as a function 

of time. T = 600 K. 
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(a) First layer. 

Figure 11.- Variation in potential 
second layers within 


Trap potential, eV 
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(b) Second layer. 

energy around a trap for first and 
the 400-site array. 


around the site is shown in figure 12 for a constant deposition rate 

= 5 X cm“2-sec*"^ and a substrate temperature T = 550 K. (Lighter shad- 

ing of some of the SOS squares represents the next adlayer of the growing film.) 

As atoms impinge on the surface, heterogeneous nucleation occurs, and growth is 
so rapid in this vicinity that captured adatoms cause a depletion zone around the 
growing cluster. Further, the vertical growth rate is also larger than normal 
because the potential energy at the trap sites is still lower on the surface of the 
growing island than on the uniform-surface sites, thus generating the early growth of 
the second layer. It is therefore likely that growth in the vicinity of such a site 
would be dominated by this defect. 

It has recently been observed that epitaxial growth is extremely sensitive to 
surface defect density and to the magnitude of the deposition rate (refs. 10 and 11). 
Growths have been achieved at very low temperatures on appropriate substrates only 
when the surface was relatively smooth and defect free (as determined by Kikuchi 
lines present in the reflection high-energy electron diffraction patterns). If the 
defect density is too high, then epitaxy is inhibited by the dominance of growth 
affected by the defects. If the deposition rate is too high, then even with low 
defect density the growth around the defects is so rapid that epitaxy is also lim- 
ited. Therefore, an understanding of the growth rate about different types of 
defects, as well as a knowledge of the overall surface defect density, would be help- 
ful in assessing the probability of epitaxy for a given system. 


Experiment 4 - Thin-Film Growth With Variable Deposition Rate and 

Substrate Temperature 

The representation of the growth for four different deposition rates over the 
temperature range 300 K to 600 K is presented in figures 13 to 16. The sequences in 
figure 13 show the dramatic difference in the growth of clusters due to the higher 
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(a) t = 0.5 sec. 



(d) t = 1 .4 sec. 



(b) t = 0.8 sec. 



(e) t = 1 .7 sec. 
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(c) t = 1 ,1 sec. 


(f) t = 2.0 sec. 


Figure 14.- Deposition at ^ ^ cm"'^-sec'“ . Lighter shading of SOS 

squares represents next adsorption layer of growing film. T = 500 K. 
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(a) t = 0,5 sec. 



(d) t = 1 .4 sec. 



(b) t = 0.8 sec. 


(e ) t = 1 . 7 sec. 
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(c ) t = 1 . 1 sec. 



(f) t = 2.0 sec. 


Figure 15.- Deposition at = 1 x 10^^ cm"^-sec'“^ • Lighter shading of SOS 

squares represents next adsorption layer of growing film. T = 600 K. 
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substrate temperature. This is simply a manifestation of the higher mobility of 
surface atoms. The growth of single clusters and their ultimate coalescence is shown 
in figure 14. This figure illustrates the development of a possible grain boundary. 
Although graphic representation of different growing grain orientations is not easily 
achieved with the SOS approach, the potential field surrounding a particular defect 
or nucleation site does provide some information about the probable growth orienta- 
tion, and this may be sufficient to determine if the coalescing grains will be in or 
near registry (i.e., there will be a low-angle grain boundary) or whether a high-angle 
grain boundary will be formed. 

Figure 15 shows the growth seguence that occurred when an advancing ledge was 
formed upon the coalescence of two large clusters. This behavior has also been 
observed by Weeks and Gilmer for calculated crystal growth from the melt (ref. 5). 

The sequence in figure 16 shows rapid growth; in only 1 second, the first layer is 
complete, the second layer is almost complete, the third layer has a large cluster, 
and there is even one adatom beginning the fourth layer. 

The roughness factor as a function of time for these growths is presented in 
figures 17 and 18, and provides an indication of the changing surface uniformity. 

The roughness factor is given by /i = where A and A are the surface areas 

of the growing film and the homogeneous substrate. The oscillating structure of the 
curves is a manifestation of the growing multilayers and the subsequent surface 
migration filling in the vacancies and the ledge and kink sites. At higher substrate 
temperatures, the more mobile adatoms fill in these sites faster; this tends to 
smooth out any structure in these curves. 


A plot of cluster density n^ and cluster size N 
presented in figure 19 for a deposition rate of 5 x 10^^ 
size (i > 1) have a total density of 


as a function of time is 
cm“^-sec"^ • Clusters of 


E 


i=2 


n. 

1 


n 
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Definite maxima in these curves are observed for all temperatures tested. As is 
apparent from figures 13 to 16, the decay in n^ is due to the coalescence of the 
clusters. This behavior has also been observed by Donohoe and Robins for Ag/NaCl 
(ref. 12), Hamilton and Logel for Ag/C and Pd/C (ref. 13), and Corbett and Boswell 
for Ag/MoS 2 (ref. 14). The maximum cluster size also decreases with increasing 
temperature, as was observed by Poppa for Bi/C and Ag/C (ref. 15). The most probable 
size of the clusters for the maxima at T = 300 K and t = 5.5 sec is approximately 
two to three atoms. 


The initial slope of the n curves in figure 19 gives the nucleation rate 
(ref. 16) 
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17,- Surface roughness factor for = 5 x 10^^ cm“^-sec"^ and several 
temperatures. Completed layers correspond to minima in curves. 
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Figure 18.- Surface roughness factor at T = 500 K for several deposition rates 
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Figure 19.- Cluster density as function of time with constant-cluster curves 

for 25, 50, 100, and 200 atoms. Decrease in n^ corresponds to growth 
coalescence. Rj = 5 x 10^^ cm"^-sec“^ . 



where Z is the Zeldovitch factor, a* is the capture number, is the cluster 

energy, and i* is the critical nucleus size. If it is assumed in this case that 
Zo* = 1 and i* = 1, then 



At T 




^o/^d 


and at 

J /R^ 
o' d 




300 K, the model gives 
-15 2 

2,93 X 10 cm -sec 

= 2.67 X 10 cm -sec 

T = 400 K, the model gives 
-16 2 

8.7 X 10 cm -sec 

-16 2 

: 7.9 X 10 cm -sec 


(r^ = 5 X 10^^ cm ^-sec 

fr. r -2 -In 

= 5 X 10 cm -sec J 

/ . 13-2 -U 

[R^ =5x10 cm -sec J 

fr. r -2 -In 

[R^ = 5 X 10 cm -sec J 


This is reasonably good agreement , and indicates that the critical nucleus is 
indeed i* = 1 . 


Experiment 5 - Annealing of Thin-Film Growth 

Annealing of a given deposition or growth proceeds by surface diffusion or mono- 
mer exchange between clusters (Ostwald ripening). As discussed previously, cluster 
peripheral motion can also occur causing an increase in the order of the growth. 

Figure 20(a) shows the initial growth condition = 5 x io^^ cm^^-sec""^ when 
T = 300 K. The temporal variation in growth rate for = 0 with the substrate 

temperature increased to 600 K and then to 700 K is presented in figures 20(b) 
through 20(e). For both temperatures, each island or cluster is driven to the more 
ordered arrangements shown and the number of smaller nuclei has noticeably decreased 
for the longer elapsed time. Figure 21 shows the effect of temperature on the 
average density of the clusters. This same sort of linear decrease in the cluster 
density was found by Donohoe and Robins for the system Au/NaCl (ref. 12). Annealing 
at deposition temperature did not seem to have a significant effect on cluster density 
even for large times, but when the temperature was increased, the low-density clusters 
ultimately broke up by monomer exchange to the larger, more stable clusters. 
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(a) 


13 -2 -1 

R_ = 5 X 10 cm -sec ; 
d 

T = 300 K; t = 0 sec. 




= 0; T = 600 K; 
= 3,0 sec. 


(d) = 0; T = 700 K; 

t = 3,0 sec. 




(c) R^ = 0; T = 600 K; 
t = 6,0 sec. 


(e) R^ = 0; T = 700 K; 
t = 6,0 sec. 


Figure 20,- Annealing of deposition. 
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Figure 21,- Annealing effects on cluster size. (n ) = 7.75 x io cm ; 

^ s T=300 

t = 6 sec; 0 = 0.3. 


CONCLUDING REMARKS 

A solid-on-solid (SOS) Monte Carlo computer simulation utilizing a potential- 
energy scaling technique has been used to model the initial stages of thin-film 
growth. The method employs Boltzmann-ordered statistics to simulate fluctuations in 
vibrational energy at each site in the 20 x 20 array. Subsequent events of 
adsorption, desorption, surface migration, incorporation, and substrate atom diffu- 
sion to the surface are considered, and their effect on the potential-energy field at 
each site is continually updated and recorded. The results of several computer 
experiments show consistency with the expected behavior of thin-film growth. Surface 
migration data taken at different substrate temperatures verified the activation 
energy as determined by an Arrhenius plot. Dispersed adatoms were observed to clus- 
ter into dimers and trimers and finally to form a single cluster in its most stable 
configuration. A point defect was designed by varying the interaction potential in 
the lateral (x,y) and vertical (z) directions, and served to illustrate the preferred 
growth that is known to occur at such defects. This procedure appears to be a prom- 
ising approach in the study of such defects. Experiments with varying deposition 
rate and substrate temperature showed the expected behavior of thin-film growth, from 
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nucleation and cluster growth to coalescence. The nucleation rate was found to be 
proportional to the square of the deposition rate, in agreement with the atomistic 
theory of nucleation. Finally, annealing experiments demonstrated the ordering of 
clusters as a function of time and temperature. 


Langley Research Center 
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NOMENCLATURE 


A 

A_ 


Ac 

s 

D 

E(i, j ) 

E3(i,j) 

Ei* 

F{e) 
f { e) 
G(e) 
g( e) 

AH 


ads 


AH 


cl 


AH 


des 


k 

N 


■"s 

p 

2d 

R 


2 

surface area of growing film, cm 

2 

surface area of homogeneous substrate, cm 

interatomic distance, cm 

. -3 

bulk concentration , cm 
surface concentration, cm 

2 —1 

diffusion coefficient, cm -sec 

total energy at site (i,j) on second and succeeding layers, eV 
total energy of surface layer at site (i,j), eV 
cluster energy, eV 

cumulative distribution function, dimensionless 
distribution function, dimensionless 

cumulative ordered distribution function, dimensionless 
ordered distribution function, dimensionless 

heat of adsorption, eV 

heat of sublimation from completed layer, eV 

heat of desorption, eV 

critical cluster size, dimensionless 
nucleation rate, cm"^-sec“^ 

Bo Itzmann ’ s cons tan t , eV/K 
coordination number, dimensionless 
cluster size, cm"^ 

random sample number, dimensionless 

. -2 
cluster density of size i, cm 

cluster density (monomers excluded), cm 

statistical probability 

activation barrier for surface diffusion, eV 
random number, dimensionless 



k 

s 

So 

T 

t 

At 

s 

U(i,j ) 
U, 


U 


m 




x,y,z 


Z 


a 

P 

Y 

6 

e 

e 

o 

^(n) 

r) 

e 

X 

^des 

c 

0 * 

T 


. . —2 -1 
deposition rate, cm sec 

r oughne s s f a c tor 

sticking coefficient, dimensionless 

sticking coefficient at zero coverage, dimensionless 
subs trate temperature , K 
time, sec 

sampling time , sec 
potential energy at site (irj)r eV 
threshold energy for incorporation, eV 
threshold energy for migration, eV 

potential energy at site (i,j) on uniform surface, eV 

coordinate system relative to film surface, cm 

Zeldovitch factor, dimensionless 

incremental energy change, nearest neighbor, eV 

incremental energy change, second-nearest neighbor, eV 

incremental energy change, third-nearest neighbor, eV 

thickness of film, cm 

random vibrational energy, eV 

maximum random vibrational energy, eV 

vibrational energy from ordered statistics, eV 

number of bonds, dimensionless 

coverage, dimensionless 

inverse energy, eV”^ 

— 2 —1 

desorption flux, cm -sec 
order of bonding 

cluster capture number, dimensionless 
stay time, sec 
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average mean time between hops 
vibrational period, sec 
single-bond energy, eV 
energy change at adsorption site, eV 
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